Ice streams in the Laurentide Ice Sheet: identification, characteristics and comparison to modern ice sheets

This paper presents a comprehensive review and synthesis of ice streams in the Laurentide Ice Sheet (LIS) based on a new mapping inventory that includes previously hypothesised ice streams and includes a concerted effort to search for others from across the entire ice sheet bed. The inventory includes 117 ice streams, which have been identified based on a variety of evidence including their bedform imprint, large-scale geomorphology/topography, till properties, and ice rafted debris in ocean sediment records. Despite uncertainty in identifying ice streams in hard bedrock areas, it is unlikely that any major ice streams have been missed. During the Last Glacial Maximum, Laurentide ice streams formed a drainage pattern that bears close resemblance to the present day velocity patterns in modern ice sheets. Large ice streams had extensive onset zones and were fed by multiple tributaries and, where ice drained through regions of high relief, the spacing of ice streams shows a degree of spatial self-organisation which has hitherto not been recognised. Topography exerted a primary control on the location of ice streams, but there were large areas along the western and southern margin of the ice sheet where the bed was composed of weaker sedimentary bedrock, and where networks of ice streams switched direction repeatedly and probably over short time scales. As the ice sheet retreated onto its low relief interior, several ice streams show no correspondence with topography or underlying geology, perhaps facilitated by localised build-up of pressurised subglacial meltwater. They differed from most other ice stream tracks in having much lower length-to-width ratios and have no modern analogues. There have been very few attempts to date the initiation and cessation of ice streams, but it is clear that ice streams switched on and off during deglaciation, rather than maintaining the same trajectory as the ice margin retreated. We provide a first order estimate of changes in ice stream activity during deglaciation and show that around 30% of the margin was drained by ice streams at the LGM (similar to that for present day Antarctic ice sheets), but this decreases to 15% and 12% at 12 cal ka BP and 10 cal ka BP, respectively. The extent to which these changes in the ice stream drainage network represent a simple and predictable readjustment to a changing mass balance driven by climate, or internal ice dynamical feedbacks unrelated to climate (or both) is largely unknown and represents a key area for future work to address

The glacial geomorphology of the western cordilleran ice sheet and Ahklun ice cap, Southern Alaska

During the late Wisconsinan, Southern Alaska was covered by two large ice masses; the western arm of the Cordilleran Ice Sheet and the Ahklun Mountains Ice Cap. Compared to the other ice sheets that existed during this period (e.g. the British-Irish, Laurentide and Fennoscandian ice sheets), little is known about the geomorphology they left behind. This limits our understanding of their flow pattern and retreat. Here we present systematic mapping of the glacial geomorphology of the two ice masses which existed in Southern Alaska. Due to spatially variable data availability, mapping was conducted using digital elevation models and satellite images of varying resolutions. Offshore, we map the glacial geomorphology using available bathymetric data. For the first time, we document >5000 subglacial lineations, recording ice flow direction. The distribution of moraines is presented, as well as features related to glacial meltwater drainage patterns (eskers and meltwater channels). Prominent troughs were also mapped on Alaska's continental shelf. This map provides the data required for a glacial inversion of these palaeo-ice masses

On the Reconstruction of Palaeo-Ice Sheets: Recent Advances and Future Challenges

Reconstructing the growth and decay of palaeo-ice sheets is critical to understanding mechanisms of global climate change and associated sea-level fluctuations in the past, present and future. The significance of palaeo-ice sheets is further underlined by the broad range of disciplines concerned with reconstructing their behaviour, many of which have undergone a rapid expansion since the 1980s. In particular, there has been a major increase in the size and qualitative diversity of empirical data used to reconstruct and date ice sheets, and major improvements in our ability to simulate their dynamics in numerical ice sheet models. These developments have made it increasingly necessary to forge interdisciplinary links between sub-disciplines and to link numerical modelling with observations and dating of proxy records. The aim of this paper is to evaluate recent developments in the methods used to reconstruct ice sheets and outline some key challenges that remain, with an emphasis on how future work might integrate terrestrial and marine evidence together with numerical modelling. Our focus is on pan-ice sheet reconstructions of the last deglaciation, but regional case studies are used to illustrate methodological achievements, challenges and opportunities. Whilst various disciplines have made important progress in our understanding of ice-sheet dynamics, it is clear that data-model integration remains under-used, and that uncertainties remain poorly quantified in both empirically-based and numerical ice-sheet reconstructions. The representation of past climate will continue to be the largest source of uncertainty for numerical modelling. As such, palaeo-observations are critical to constrain and validate modelling. State-of-the-art numerical models will continue to improve both in model resolution and in the breadth of inclusion of relevant processes, thereby enabling more accurate and more direct comparison with the increasing range of palaeo-observations. Thus, the capability is developing to use all relevant palaeo-records to more strongly constrain deglacial (and to a lesser extent pre-LGM) ice sheet evolution. In working towards that goal, the accurate representation of uncertainties is required for both constraint data and model outputs. Close cooperation between modelling and data-gathering communities is essential to ensure this capability is realised and continues to progress

Exploring the ingredients required to successfully model the placement, generation, and evolution of ice streams in the British-Irish Ice Sheet

Ice stream evolution is a major uncertainty in projections of the future of the Greenland and Antarctic Ice sheets. Accurate simulation of ice stream evolution requires an understanding of a number of “ingredients” that control the location and behaviour of ice stream flow. Here, we test the influence of geothermal heat flux, grid resolution, and bed hydrology on simulated ice streaming. The palaeo-record provides snapshots of ice stream evolution, with a particularly well constrained ice sheet being the British-Irish Ice Sheet (BIIS). We implement a new basal sliding scheme coupled with thermo-mechanics into the BISICLES ice sheet model, to simulate the evolution of the BIIS ice streams. We find that the simulated location and spacing of ice streams matches well with the empirical reconstructions of ice stream flow in terms of position and direction when simple bed hydrology is included. We show that the new basal sliding scheme allows the accurate simulation for the majority of BIIS ice streams. The extensive empirical record of the BIIS has allowed the testing of model inputs, and has helped demonstrate the skill of the ice sheet model in simulating the evolution of the location, spacing, and migration of ice streams through millennia. Simulated ice streams also prompt new empirical mapping of features indicative of streaming in the North Channel region. Ice sheet model development has allowed accurate simulation of the palaeo record, and allows for improved modelling of future ice stream behaviour

An updated radiocarbon-based ice margin chronology for the last deglaciation of the North American Ice Sheet Complex

The North American Ice Sheet Complex (NAISC; consisting of the Laurentide, Cordilleran and Innuitian ice sheets) was the largest ice mass to repeatedly grow and decay in the Northern Hemisphere during the Quaternary. Understanding its pattern of retreat following the Last Glacial Maximum is critical for studying many facets of the Late Quaternary, including ice sheet behaviour, the evolution of Holocene landscapes, sea level, atmospheric circulation, and the peopling of the Americas. Currently, the most up-to-date and authoritative margin chronology for the entire ice sheet complex is featured in two publications (Geological Survey of Canada Open File 1574 [Dyke et al., 2003]; ‘Quaternary Glaciations – Extent and Chronology, Part II’ [Dyke, 2004]). These often-cited datasets track ice margin recession in 36 time slices spanning 18 ka to 1 ka (all ages in uncalibrated radiocarbon years) using a combination of geomorphology, stratigraphy and radiocarbon dating. However, by virtue of being over 15 years old, the ice margin chronology requires updating to reflect new work and important revisions. This paper updates the aforementioned 36 ice margin maps to reflect new data from regional studies. We also update the original radiocarbon dataset from the 2003/2004 papers with 1541 new ages to reflect work up to and including 2018. A major revision is made to the 18 ka ice margin, where Banks and Eglinton islands (once considered to be glacial refugia) are now shown to be fully glaciated. Our updated 18 ka ice sheet increased in areal extent from 17.81 to 18.37 million km2, which is an increase of 3.1% in spatial coverage of the NAISC at that time. Elsewhere, we also summarize, region-by-region, significant changes to the deglaciation sequence. This paper integrates new information provided by regional experts and radiocarbon data into the deglaciation sequence while maintaining consistency with the original ice margin positions of Dyke et al. (2003) and Dyke (2004) where new information is lacking; this is a pragmatic solution to satisfy the needs of a Quaternary research community that requires up-to-date knowledge of the pattern of ice margin recession of what was once the world’s largest ice mass. The 36 updated isochrones are available in PDF and shapefile format, together with a spreadsheet of the expanded radiocarbon dataset (n = 5195 ages) and estimates of uncertainty for each interval

Data for: Reconciling records of ice streaming and ice margin retreat to produce a palaeogeographic reconstruction of the deglaciation of the Laurentide Ice Sheet

The attached Excel spreadsheet provides additional information for individual ice streams covered in the manuscript 'Reconciling records of ice streaming and ice margin retreat to produce a palaeogeographic reconstruction of the deglaciation of the Laurentide Ice Sheet

Reconciling records of ice streaming and ice margin retreat to produce a palaeogeographic reconstruction of the deglaciation of the Laurentide Ice Sheet

This paper reconstructs the deglaciation of the Laurentide Ice Sheet (LIS; including the Innuitian Ice Sheet) from the Last Glacial Maximum (LGM), with a particular focus on the spatial and temporal variations in ice streaming and the associated changes in flow patterns and ice divides. We build on a recent inventory of Laurentide ice streams and use an existing ice margin chronology to produce the first detailed transient reconstruction of the ice stream drainage network in the LIS, which we depict in a series of palaeogeographic maps. Results show that the drainage network at the LGM was similar to modern-day Antarctica. The majority of the ice streams were marine terminating and topographically-controlled and many of these continued to function late into the deglaciation, until the ice sheet lost its marine margin. Ice streams with a terrestrial ice margin in the west and south were more transient and ice flow directions changed with the build-up, peak-phase and collapse of the Cordilleran-Laurentide ice saddle. The south-eastern marine margin in Atlantic Canada started to retreat relatively early and some of the ice streams in this region switched off at or shortly after the LGM. In contrast, the ice streams draining towards the north-western and north-eastern marine margins in the Beaufort Sea and in Baffin Bay appear to have remained stable throughout most of the Late Glacial, and some of them continued to function until after the Younger Dryas (YD). The YD influenced the dynamics of the deglaciation, but there remains uncertainty about the response of the ice sheet in several sectors. We tentatively ascribe the switching-on of some major ice streams during this period (e.g. M'Clintock Channel Ice Stream at the north-west margin), but for other large ice streams whose timing partially overlaps with the YD, the drivers are less clear and ice-dynamical processes, rather than effects of climate and surface mass balance are viewed as more likely drivers. Retreat rates markedly increased after the YD and the ice sheet became limited to the Canadian Shield. This hard-bed substrate brought a change in the character of ice streaming, which became less frequent but generated much broader terrestrial ice streams. The final collapse of the ice sheet saw a series of small ephemeral ice streams that resulted from the rapidly changing ice sheet geometry in and around Hudson Bay. Our reconstruction indicates that the LIS underwent a transition from a topographically-controlled ice drainage network at the LGM to an ice drainage network characterised by less frequent, broad ice streams during the later stages of deglaciation. These deglacial ice streams are mostly interpreted as a reaction to localised ice-dynamical forcing (flotation and calving of the ice front in glacial lakes and transgressing sea; basal de-coupling due to large amount of meltwater reaching the bed, debuttressing due to rapid changes in ice sheet geometry) rather than as conveyors of excess mass from the accumulation area of the ice sheet. At an ice sheet scale, the ice stream drainage network became less widespread and less efficient with the decreasing size of the deglaciating ice sheet, the final elimination of which was mostly driven by surface melt

Ribbed bedforms on palaeo-ice stream beds resemble regular patterns of basal shear stress (‘traction ribs’) inferred from modern ice streams

Rapidly-flowing ice streams are an important mechanism through which ice sheets lose mass, and much work has been focussed on elucidating the processes that increase or decrease their velocity. Recent work using standard inverse methods has inferred previously-unrecognised regular patterns of high basal shear stress (‘sticky spots’ >200 kPa) beneath a number of ice streams in Antarctica and Greenland, termed ‘traction ribs’. They appear at a scale intermediate between smaller ribbed moraines and much larger mega-ribs observed on palaeo-ice sheet beds, but it is unclear whether they have a topographic expression at the bed. Here, we report observations of rib-like bedforms from DEMs along palaeo-ice stream beds in western Canada that resemble both the pattern and dimensions of traction ribs. Their identification suggests that traction ribs may have a topographic expression that lies between, and partly overlaps with, ribbed moraines and much larger mega-ribs. These intermediate-sized bedforms support the notion of a ribbed bedform continuum. Their formation remains conjectural, but our observations from palaeo-ice streams, coupled with those from modern ice masses, suggest they are related to wave-like instabilities occurring in the coupled flow of ice and till and modulated by subglacial meltwater drainage. Their form and pattern may also involve glaciotectonism of subglacial sediments

Ice stream activity scaled to ice sheet volume during Laurentide Ice Sheet deglaciation

The contribution of the Greenland and West Antarctic ice sheets to sea level has increased in recent decades, largely owing to the thinning and retreat of outlet glaciers and ice streams1, 2, 3, 4. This dynamic loss is a serious concern, with some modelling studies suggesting that the collapse of a major ice sheet could be imminent5, 6 or potentially underway7 in West Antarctica, but others predicting a more limited response8. A major problem is that observations used to initialize and calibrate models typically span only a few decades, and, at the ice-sheet scale, it is unclear how the entire drainage network of ice streams evolves over longer timescales. This represents one of the largest sources of uncertainty when predicting the contributions of ice sheets to sea-level rise8, 9, 10. A key question is whether ice streams might increase and sustain rates of mass loss over centuries or millennia, beyond those expected for a given ocean–climate forcing5, 6, 7, 8, 9, 10. Here we reconstruct the activity of 117 ice streams that operated at various times during deglaciation of the Laurentide Ice Sheet (from about 22,000 to 7,000 years ago) and show that as they activated and deactivated in different locations, their overall number decreased, they occupied a progressively smaller percentage of the ice sheet perimeter and their total discharge decreased. The underlying geology and topography clearly influenced ice stream activity, but—at the ice-sheet scale—their drainage network adjusted and was linked to changes in ice sheet volume. It is unclear whether these findings can be directly translated to modern ice sheets. However, contrary to the view that sees ice streams as unstable entities that can accelerate ice-sheet deglaciation, we conclude that ice streams exerted progressively less influence on ice sheet mass balance during the retreat of the Laurentide Ice Sheet

Ice streams in the Laurentide Ice Sheet: a new mapping inventory

Rapidly flowing ice streams dominate the drainage of continental ice sheets and are a key component of their mass balance. Due to their potential impact on sea level, their activity in the Antarctic and Greenland Ice Sheets has undergone detailed scrutiny in recent decades. However, these observations only cover a fraction of their ‘life-span’ and the subglacial processes that facilitate their rapid flow are very difficult to observe. To circumvent these problems, numerous workers have highlighted the potential of investigating palaeo-ice streams tracks, preserved in the landform and sedimentary record of former ice sheets. As such, it is becoming increasingly important to know where and when palaeo-ice streams operated. In this paper, we present a new map of ice streams in the North American Laurentide Ice Sheet (LIS; including the Innuitian Ice Sheet), which was the largest of the ephemeral Pleistocene ice sheets and where numerous ice streams have been identified. We compile previously published evidence of ice stream activity and complement it with new mapping to generate the most complete and consistent mapping inventory to date. The map depicts close to three times as many ice streams (117 in total) compared to previous inventories, and categorises them according to the evidence they left behind, with some locations more speculative than others. The map considerably refines our understanding of LIS dynamics, but there is a clear requirement for improved dating of ice stream activity